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The natural explanation of the excitation spectrum of the spin- 
1 antiferromagnetic Heisenberg chain is given from the viewpoint 
of the spin- 1/2 ferromagnetic-antiferromagnetic alternating Heisen- 
berg chain. The energy spectrum of the latter is calculated with 
fixed momentum k by numerical diagonalization of finite size sys- 
tems. It consists of a branch of propagating triplet pair (triplet 
wave) and the continuum of multiple triplet waves for weak fer- 
romagnetic coupling. As the ferromagnetic coupling increases, the 
triplet wave branch is absorbed in the continuum for small k : repro- 
ducing the characteristics of the spin-1 antiferromagnetic Heisen- 
berg chain. 
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Recently, the spin-1/2 ferromagnetic-antiferromagnetic alternat- 
ing Heisenberg chain ffT|— |T0|1 (ALHC) has been studied extensively to 
clarify the nature of the Haldane state of the spin-1 antiferromag- 
netic Heisenberg chain (AFHC) [|TTHZ8j|. In the present work, we 
further investigate the excitation spectrum of ALHC to understand 
the physical nature of the excitation spectrum of the spin-1 AFHC. 

The Hamiltonian of the spin-1/2 ALHC is given as follows: 

N N 
H = 2 J ^2 ^2x^21+1 — 2 J S2i-\S2ii (1) 

where Si is the spin 1/2 operator on the i-th site and the number 
of the lattice sites is 2N. We assume J, J' > and denote the 
ratio J'/ J by (3. For (3 — » oo, the spins S2i-i and S^i form a local 
triplet and this model reduces to the spin-1 AFHC whose excitation 
spectrum is speculated as follows: 

The lowest excitation has the wave number k = tt (zone bound- 
ary). The excitations near k ~ 7r form a single magnon branch 
which is not only separated from the ground state by the Haldane 
gap but also from the higher energy continuum by a finite gap. On 
the other hand, near k ~ 0, the lowest excitation is the lower edge 
of the continuum of the two magnon scattering states. At k = 0, 
the excitation energy is nearly twice the Haldane gap as expected 
from the above picture. The schematic excitation spectrum is sum- 
marized in Fig. [T](a). This argument is not only based on the 
numerical results [|I3}-|I7j| but also is supported by the original Hal- 
dane conjecture based on the nonlinear cr-model, because the only 
single particle excitations of the nonlinear cr-model are spin triplet 
bosons [PTTHT^I whose mass (Haldane gap) is much smaller than its 
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band width. The excitation spectrum of the spin-1 AFHC is also 
studied by the neutron scattering experiments for real materials 



For small (3, the excitation spectrum of the spin- 1/2 ALHC may 
be understood using the lowest order perturbation theory in (3. As 
discussed by the author in the previous work , the lowest excited 
state is the propagating wave of a triplet pair (triplet wave; hereafter 
abbriviated as TW) which breaks one of the local singlets on the 
antiferromagnetic bonds. The excitation energy E{k) of this state 
with momentum k is given by 

£(fc)~2J[l + ^cos(fc)], (2) 

where —tt < k < tt and the length of the unit cell(= twice lattice 
constant) is set equal to unity. This type of excitation has been in- 
troduced by Knabe [[n| for the AKLT model [|?Uj • Fath and Solyom 
|2"If1 have shown that this excitation is essentially equivalent to the 
domain wall defect in the long range string order [|^4|^5|1 introduced 
by Arovas, Auerbach and Haldane [|24|J25| . The same is true for the 



present case. The domain wall character of the elementary excita- 
tion in the Haldane state is also verified numerically for the spin-1 
AFHC by Sakai and Takahashi JTB 



The excitations with higher energies are made up of multiple 
TW excitations which form a continuum above the TW branch. 
The energy of the excitation with momentum k in the two TW 
continuum is given by 

E(k\ q) ~ 4 J[l + ^{cos(g) + cos(/c - q)}} (3) 
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with arbitrary value of q. Thus, the upper and lower boundaries of 
the continuum at momentum k are given as follows, 

4J[1 - |cos(fc/2)] < E{k;q) < 4J[1 + | cos(fc/2)]. (4) 

The schematic energy spectrum is shown in Fig. [3(b). At first sight, 
this spectrum appears to be quite different from Fig. [3(a). In the 
following, we examine how the spectrum of Fig. 13(b) changes into 
that of Fig. 0(a) with the increase of (3. 

Within the lowest order perturbation theory in (3, we may specu- 
late as follows: The gap between the TW branch and the continuum 
is minimum (= 2J(1 — \(3)) at k = and decreases with the increase 
of (3. This gap would collapse around (3 = (3 C ~ | and for (3 > (3 C the 
TW branch would be absorbed in the two TW continuum for small 
k. Then, the lowest excitation at k = would be replaced by the 
lower edge of the continuum as observed in the spin-1 antiferromag- 
netic chain. We verify this senario by the numerical diagonalization 
of finite size systems. 

We have diagonalized the finite-size Hamiltonian ([3) for 2N = 
8, 12, 16, 20 and 24 with the periodic boundary condition. We have 
modified the program package TITPACK version 2 supplied by 
Nishimori to diagonalize the Hamiltonian in the subspace with fixed 
momentum k. The excitation energies E(k) of the low lying ex- 
cited states with momentum k are obtained. The excitation spec- 
tra are shown in Fig. 0(a)- (d) for 2N = 24. For small (3 (Fig. 
0(a)), the low energy excitation spectrum is qualitatively similar to 
Fig. [3(b). Namely, it consists of a TW branch separated from the 
higher energy states by an energy gap of the order of J. This gap 
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decreases with the increase of (3 and finally collapses at a critical 
value (3 = (3 C ~ 1.5 (Fig.|2|(b)). We estimate the value of (3 C more 
precisely afterwards. 

In order to get more insight about what happens for large /3, 
we classify the states by the total spin 5 tot which is also a good 
quantum number in the isotropic case. As shown in Fig. g, the 
lowest excited states with non-zero momentum have always S tot = 1 
while those with k = have S' tot = 1 for small /3 and S tot = 2 for 
large (3. 

In Figs. |3](a-d), we plot the system size dependence of 
E (0,1,N), E (0,2,N) and 2E (ir, 1, N) (twice the Haldane gap) 
for 2N = 8, 12, 16, 20 and 24 where E (k, S to \ N) denotes the low- 
est excitation energy of the 2N membered chain with momentum 
k and total spin S tot . These data are extrapolated to N — > oo by 
means of the Shanks' transform (j29j. The quantities Eq(0,2,N) 
and 2Eq(7t, 1, N) converge to the almost same value for all f3. This 
means that the lowest excited state with S = 2 is the scattering 
states of two TW's with momenta tt. 

For f3 > (3 C , Eq(0, 1,N) also converges to the same value (Figs. 
|3|(c,d)) while it converges to a different lower value for (3 < (3 C (Figs. 
|3|(a,b)). Therefore the lowest state with 5 tot = 1 at k = is the 
scattering state two TW's or a single TW state according as (3 > (3 C 
or < (3 C . This implies that the long wave length part of the TW 
branch is already absorbed in the continuum for (3 > (3 C including 
the limit of spin-1 AHFC. 

As discussed by White and Huse [jT5j] for spin-1 AFHC, two 
triplet states with identical set of quantum numbers cannot form a 
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totally triplet state due to symmetry. This means that two TW's 
with k = 7r cannot form a four spin triplet state with k = 0. There- 
fore, the two TW state with S tot = 1 at k = must be the scattering 
state of two TW's with different momenta slightly deviated from w. 
This explains why the state with S tot = 1 has higher energy than 
that with S tot = 2 at k = for (3 > (3 C in the finite size system. 

Making use of this property, we may determine the critical value 
(3 C more precisely by extrapolating the crossing point (3 C (N) of 
E (0,1,N) and E (0,2,N) to N -> oo. The system size depen- 
dence of (3 C (N) is plotted in Fig.f|. Assuming the asympototic form 
P C {N) ~ p c + aA^, we find (3 C ~ 1.62. 

To summarize, the excitation spectrum of the spin-1/2 alternat- 
ing Heisenberg chain is calculated by the lowest order perturbation 
in P and the numerical diagnalization of finite size systems. The na- 
ture of the lowest excited state at k = change from the single TW 
state to the two TW state around (3 = j3 c ~ 1.62. The characteristic 
features of the excitation spectrum of the spin-1 antiferromagnetic 
Heisenberg chain is naturally understood as the large (3 limit of the 
spin-1/2 alternating Heisenberg chain. 

The author is grateful to H. Nishimori for providing TITPACK 
version 2. This work is supported by the Grant-in- Aid for Scien- 
tific Research on Priority Areas " Computational Physics as a New 
Frontier in Condensed Matter Research" from the Ministry of Edu- 
cation, Science and Culture. The numerical calculation is performed 
by HITAC S3800/480 at the Computer Center of the University of 
Tokyo and HITAC S820/15 at the Information Processing Center 
of Saitama University. 



6 



References 

fi [1]K. Hida: Phys. Rev. B45 (1992) 2207. 

fi [2]K. Hida: Phys. Rev. B45 (1992) 8268; see also 

the comment by K. Okamoto and D. Nishino: Phys. Rev. B48 
(1993) 13149. 

fi [3]K. Hida and S. Takada: J. Phys. Soc. Jpn. 61 (1992)1879. 

fi [4]K. Hida: J. Phys. Soc. Jpn. 62 (1993) 439; ibid. 1463; ibid. 
1466; ibid. 3357. 

fi [5]M. Kohmoto and H. Tasaki: Phys. Rev. B46 (1992) 3486. 

fi [6]M. Yamanaka, Y. Hatsugai and M. Kohmoto: Phys. Rev. 
B48 (1993) 9555. 

fi [7]S. Takada: J. Phys. Soc. Jpn. 61 (1992) 428. 

fi[8]N. Hatano and M.Suzuki: J. Phys. Soc. Jpn. 62 (1993) 847. 

fi [9]K. Okamoto, D. Nishino and Y. Saika: J. Phys. Soc. Jpn. 
62 (1993) 2587. 

fi [10]K. Okamoto and D. Nishino: preprint (1994). 

7 



fi [11JF.D.M. Haldane: Phys. Lett. 93A (1983) 464; Phys. Rev. 
Lett. 50 (1983) 1153. 

fi [12]I. Affleck: J. Phys. Condens. Matter. 1 (1989) 3047 and 
references therein. 

fi[13]I. Affleck and R. A. Weston: Phys. Rev. B45 (1992) 4667. 

fi [14] E. S. S0rensen and I. Affleck: Preprint 

(cond-mat.9311021)(1993). 

fi [15]S.R. White and DA. Huse: Phys. Rev. B48 (1993) 3844. 

fi [16] O. Golinelli, Th. Jolicoeur and D. R. Lacaze: Phys. Rev. 
B46 (1992) 10854. 

fi [17] M. Takahashi: Phys. Rev. Lett. 62 (1989) 2313; Phys. 
Rev. B48 (1993) 311. 

fi [18]T. Sakai and M. Takahashi: J. Phys. Soc. Jpn. 63 (1994) 
755. 

fi [19]S. Knabe: J. Stat. Phys. 52 (1988) 627. 

fi [20] I. Affleck, T. Kennedy, E. Lieb and H. Tasaki: Phys. Rev. 

8 



Lett. 59 (1987) 799; Commun. Math. Phys. 115 (1988) 477. 

fi [21]G. Fath and J. Solyom: J. Phys. Condens, Matter. 5 
(1993) 8983. 

fi [22] M. den Nijs and K. Rommelse: Phys. Rev. B40 4709 
(1989). 

fi [23]H. Tasaki: Phys. Rev. Lett. 66 798 (1991). 

fi [24]D. P. Arovas, A. Auerbach and F. D. M. Haldane: Phys. 
Rev. Lett. 60 (1988) 531. 

fi [25] M. Kaburagi, I. Harada and T. Tonegawa: J. Phys. Soc. 
Jpn. 62 (1993) 1848. 

fi [26] S. Ma, C. Broholm, D. H. Reich, B. J. Sternlieb and R. 
W. Erwin: Phys. Rev. Lett. 69 (1992) 3571. 

fi [27]J.P. Renard, M. 

Verdaguer, L.P. Regnault, W.A.C. Erkelens, J. Rossat-Mignod 
and W.G. Stirling: Europhys. Lett. 3 (1987) 945. 

fi [28] J. P. Renard, M. Verdaguer, L.P. 



9 



Regnault, W.A.C. Erkelens, J. Rossat-Mignod, J. Ribas, W.G. 
Stirling and C. Vettier: J. Appl. Phys. 63 (1988)3538. 

fi [29]D. Shanks: J. Math. Phys. 34 (1955) 1. 



10 



FIGURES 

Fig. 1. Schematic excitation spectrum of the spin-1 antiferro- 
magnetic Heisenberg chain (a) and the spin- 1/2 alternating Heisen- 
berg chain with small (3 (b). 

Fig. 2. Excitation spectrum of the alternating Heisenberg chain 
with 2N = 24 with (3 = 0.5 (a), 1.5 (b), 2.0 (c) and 4.0(d). The 
values of total spin S tot are 0(o), 1 (•), 2 (□) and 3 (O). Energy 
unit is J. 

Fig. 3. System size dependence of the lowest excitation energy 
with k = and S tot = 1(»), k = and S tot = 2(D) and twice 
the Haldane gap with k = n and S tot = 1(A) for (3 = 0.5 (a), 1.0 
(b), 1.5 (c) and 2.0(d). The values extrapolated to TV — > oo by the 
Shanks' transform are also shown. Energy unit is J. 

Fig. 4. System size dependence of the crossing point (3 C (N) of 
E (0, 1, N) and E (0, 2, N) (•). The solid line is the least square fit 
to the asympototic form (3 C (N) ~ (3 C + aN~ v with (3 C ~ 1.62. 
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